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Abstract 



The 3d Ising model on a regular cubic lattice can be expressed in terms of an SU (2) 2d 
fermionic model with a Z2-fiuxes. We modify the model such that it is defined on the dual 
to a body centered cubic lattice. The advantage of this lattice is that 2d embedded surfaces 
have no selfintersections, thus partially avoiding the Sign-factor problem associated with 
the 2d fermionoc models related to the 3d Ising model. Rather than solving the full SU{2) 
fermionic theory on this lattice we study the simpler model of scalar fermions and find the 
spectrum of excitations. The model has no mass gap. We reformulate the model using the 
i?-formalism and a new interesting structure appears due to the necessity of introducing 

(3) 

a three-paticle matrix -R)^-^. It encodes the essential character of the Sign-factor. We 
analyse the integrability properties of this class of models. 
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1 Introduction 



After a conjecture by A.Polyakov [|^ that the three-dimensional Ising Model (3DIM) 
can be represented as a fermionic string theory in a three dimensional Euclidean space 
a number of works was carried out during the eighties p|, H, ||, |], |, ^, |10| trying 
to substantiate the conjecture. The main problem is the so called Sign-factor, which is 
the analogy of the Kac-Ward factor [11| for the Ising Model formulated on a 2D regular 
lattice. This Sign-factor ensures the cancellation of self-intersecting 2D surfaces mapped 
onto the target space 3D regular lattice. The cancellation is necessary in order to have 
correct Boltzmann weights in the string representation of the 3D statistical sum as an 
ensemble of the random 2D surfaces 0, ^. It resembles the Pauli principle for 
Fermi particles, which prohibits two of them to be in the same state (in the same space- 
time point in this case). Therefore it is expected that the Sign-factor is connected with 
some fermionic structure on the 2d world-sheet of the strings. 

In PI a fermionic model with the Sign-factor was defined on the 2D regular lattice. 



and its continuum limit was investigated in ||^, |T2|. The model appeared to be a model 
of SO (3) fermions hopping in the staggered field of Z2-fluxes on a random Manhattan 
lattice. 

A similar model of fermions hopping in a staggered f/(l)-gauge field background on a 
regular Manhattan lattice was analysed in [|T^. 

In ||15[ a model of spins with Z2 global symmetry was considered on the so called dual 
Body Centered Cubic (BCC) lattice (see Fig.l). The BCC lattice consists of two simple 
cubic sublattices arranged in such a way that the sites of one sublattice are positioned in 
the middle points of the cubes of the second sublattice (see Fig. la). 

The dual to the BCC (DBCC) lattice is represented in the Fig. lb. As one can see, 
the two dimensional faces of the dual BCC lattice are: 

i) hexagons, which are dual to the links connecting the neighbour vertices of different 
sublattices (dotted hues in the Fig. lb) and 

ii) squares, which are orthogonal to the links of the same sublattice (bold lines in the 
Fig. lb). 

This lattice is interesting in a sense that the relevant surfaces here, i..e. those which 
at most occupy the faces of the DBCC lattice once, have no self-intersections at all. The 
analogy for contours on a 2D lattice is the honeycomb (or its dual triangular) lattice, 
where curves have no self-intersection points. 

We follow the same reasoning as in in the case of a regular lattice, but now on a 
DBCC lattice. We then obtain a fermionic hopping model in a staggered SU (2)-gauge 
field, defined on the Manhattan lattice of the type shown in Fig. 3. On this lattice we do 
not have the magnetic fluxes which represented the ends of the self-intersection lines on 
the cubic lattices. However, we will see the emergence of a new structure of the action, 
which encodes the essential character of the Sign-factor and which is interesting from the 
point of view of integrability of the model. 



1 



In this article we define and investigate a model of fermions hopping in a staggered 
?7(l)-gauge field background on the ML defined in Fig.3. We find the spectrum of the 
model and analyse its continuum limit. We introduce a three particle scattering R^- 
matrix, which, together with the ordinary two particle i?^-matrix constitutes the Mon- 
odromy matrix of the model by appearing alternating in a product. We analyse the inte- 
grability property of this type of model and write down the corresponding Yang-Baxter 
equations (YBE). 

2 The Model 

We consider now a simple 2d fermionic system hopping on the ML (Fig.3) which 
was constructed on particular surfaces on the DBCC lattice as mentioned above. This 
simplification will allow us to analyse the continuum limit of the model and investigate 
integrability structure. 

The ML structure defines the hopping of fermions only along the arrows of the lattice 
(see Fig. 3) and the Hamiltonian is non hermitian: 

-f^ = ^ tH,rnC^in)UH,rnc{m) + ^c^{n)c{n). (1) 

<n,m> n 

Here Ufi^rn are the group elements of the external f/(l) field defined according to 
demand that the magnetic fluxes around the all plaquettes with circulating arrows will be 
= 27ip/q (p and q are mutually prime integer numbers). The hopping parameters t^^^ 
are chosen to be periodic, reflecting the translational symmetry of the lattice by vectors 
2ix, 2(1 + \/3)iy, where and iy are unit vectors in coordinates direction. It allows us to 
distinguish 10 types of particles (Fig.l): 

In order to diagonalize the Hamiltonian (0) we pass to the Bloch wave basis: 



cAk) = , ^ e 



'^<y) - ~irj ' - U, ...,iV^.(y) - i. 

iVa;(j,) 

{Nx,2Ny are the numbers of the hexagons in the corresponding directions), and the 
Hamiltonian becomes: 

H{U„^) = J24ik)H,,ik,^)c,in), (2) 
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where H{kx, ky) is a 10 x 10 matrix: 



^ 1 ti2e»P ti4e'9 ^ 

taie-'P 1 ^ ^ 

120 e ^ 

t32e-*9 1 t34e-*P 

t43e-»P 1 t45e*<3 

, 1 t56e"3 

1 t67e-»P i69e*9 

, iP tree-ip 10 

tsTe-'? 1 tsge^P 

tgse^P 1 ^ ^ 

\ ^ ^ ^ ^ I I 

\ toie 2 lose ' 



We have defined fc^; = p, -y/S/Cy = q, p + q = Q and p — q = P. 

The partition function for the excitations of the states with energy E is the functional 
integral over the Grassmann variables ipiln) , ijji{n) corresponding to the fermions. 

Z{E) = tre-^'' " / n d^i{n)diJiin)e-''^^''''''^+^^^ . (3) 
Since this integral is Gaussian with respect to ipiln) , tpiln) fields it follows that 

Z{E)= II [detH{k^,ky)-E]. (4) 

The equation detiJ(A;^, ky) = gives the critical line for the low lying excitations. 

The investigation of the excitations is conveniently done by passing from the La- 
grangian (action) formalism in the (1 + 1) Euclidean space-time (formula ^ to the Hamil- 
tonian and consider the representation of the partition function via the Transfer Matrix 
as follows 

Z{0) = tr[Tfy. (5) 

By analysing the periodic structure of the ML under consideration, the complete 
Transfer Matrix can be written as the product of four transfer matrices Tj, i = 1,2,3,4 
each of which mediates the evolution between the states \t + i — l) and \t + i), respectively, 
(Fig. 3): 

\t + i) =T,\t + i-l) (6) 
and the complete Transfer Matrix is equal to 

T = nr^nn. (7) 
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The technique of passing from the Lagrangian to the Hamiltonian for models on ML 
was developed in [0, based on the concept of coherent states 0]. It is possible to intro- 
duce two types of coherent states which are the eigenstates of the creation or annihilation 



operators, respectively. As shown in [0 for a ML, we need to consider two types of 
coherent sates alternating in a chain. One state differs from other by particle-hole trans- 
formation. Consequently we will need a particle and hole ordering prescription in order 
to complete the definition of the Transfer Matrices. Let us consider particle and hole 
\%b) states as follows 



1^2.) 
I^2i+l) 
(V'2i|^2i) 

with the properties 



(4i+l - ^2i+l)|0), (^2i+l| = (0|(c2i+l + i)2i+l) 



d^p2i+ld^p2^+S2i+l){i^2^+l\e^''+''^'' 



1, 



(9) 



As it is clear from Fig. 3, the Transfer Matrices T2 and T4 are permuting the disposition 
of particle and hole type of coherent states at the odd and even sites, and because of this 
the definition of the coherent state of the whole chain at the Euclidean times \t+2) , |t+3) 
must be interchanged. 

Following it is easy to see that the Transfer Matrices are of the form : 

Ti=:exp-Hr., (10) 
Hi = ^[-ti2C+(2i - l)c2(2i) + t344(2i)ci(2i - 1) - 



(tl4 + 1)C+(2Z - l)ci{2l - 1) + (t32 + 1)4(22)C2(22)], (11) 
^[-^434(2^ + 1)C2(2Z) + t45i534(2^ " ^M^l) + ^764(2^ + l)Cl{2t) - 



h^h-jct{2i - l)ci(22) - (1 - t56t45)4(2^)ci(2z - 1) + 

(1-^75^53)4(2^02(22-1)], 

^[-t674(2^)c2(22 - 1) + t894(2^ - l)ci(20 + 



(12) 



(t69 + l)4(20ci(22) + (% + 1)4(2^ - 1)C2(22 - 1)], (13) 
^[-^sct{2i)c2{2i + 1) + t9oto84(2«)c2(22 - 1) + t2ict{2i)ci{2i + 1) - 



- t2otoi4(2^)ci(2« - 1) + (1 - t9oto,i)4(20ci(2^ - 1) + 

+ (l-t20t08)4(20c2(2z-l)]. 



(14) 
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where under the simbol :: we mean ordinary normal ordering of the fermionic operators 
at the even sites of the chain and anti-normal ordering for the odd sites. 

Evaluating now the product T1T2T3T4 for the whole Transfer Matrix T (|^) in the basis 
of coherent states and puting it into the formula (H), one can obtain the expression (H) 
for the partition function Z. 

Making a Fourier transformation of the Hamiltonians ([T0|) and expressing the product 
(0) of the Transfer matrices by use of the generators of the sli algebra in a Schwinger 
form 

T = De"-^ = D JJ e"^^ (15) 

V 

Hp = eSp + n{npi+np2), (16) 

5*1 = C1C2-C2C1, S2 = C1C2 + C2C1, (17) 

S3 = cfci - 4c2, So = I, (18) 

one can, after some algebra, obtain the spectrum of the excitations. In order to do that, 
and for simplicity, from now on we will consider two particular types of parameterization 
of the hopping amplitudes: 



i) 



tl2 — tlA — ^34 — ^32 — ^67 — ^69 — ^87 — ^76 — ^89 — ^98 — ^21 — ^Ai — ^ 



1 

1 



^53 — ^56 — ^45 — ^75 — ^08 — ^90 — ^20 — ^01 — (19) 



and 
ii) 



^iO ^Oj ^1 (0C\\ 

tij=ti, z,j7^5,0, ti/t = x. 

The first parameterization is chosen by simplicity arguments while the second is demanded 
by the rotation invariance. 

First we will analyse the case i). The eigenvalues of the Hamiltonian ([ISD is found 
to he E = ±e(p) with the following dispersion relation, chemical potential fi and the 
parameter D 

cosh e(p) = 2 sin^ f ((1 + 4 cos^ $) + 4 cos - 2 cos $) + 2 + 1 /2t-^ + 2 cos ^t'^ 
+2 cos 2p (2 sin^ |((2 cos $ - 1)^^ - 2 cos + 1) - 2 cos $ - t~^) 
+ cos4p(l -4sin2 ft^) , 



fi =0, 
D =t^^^ 

The condition when the massless excitations appear is defined by the equation 

2 sin^ f ((1 - 8 sin^ f cos $)t8 + (8 cos $ - 4)^^ + 4 sin' f t^) + i = 
= 2t4(2cos$-l)+4sin2ft2. 



(21) 



(22) 
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The limit t oo, $ ^ on the critical line (P^ ) becomes 

sm|t' = l, (23) 

and the dispersion relation simplifies to 

cosh e = 12(1 - cos 2p) + cos 4p. (24) 

We see that at the point p=0 there is no gap in the spectrum £(0) = 0, and the 
excitations energy near that point is a linear function of the momentum: 

e{p) = ±AV2p. (25) 

In the case ii) the model under discussion here is connected with the Chalker-Coddington 
phenomenological model [0, |13| for the edge excitations in Hall effect responsible for the 
plateau-plateau transitions. 

The equation of the spectrum in this case is the following 

coshe(p) = 1 + (cos4p - l)(t"^ + 2(cos<l> - x^)t) + 

2(cos 2p - 1) {xt^ cos $(2x2 cos $ - 1 - x^)+ (26) 
—2t'^x~^ cos <l>(cos $ — x^) + a;(a;2 — cos $) — cos $ — t^^x~^) , 

and the equation of the critical line takes the form 

{x^ + 1- 2x^ cos^)x^ (2ti° cos $ — 2x cos + x'^t^/2) 

+ {x^ - COS $) {At^x cos $ - 2fx^ -t^ + 2t^x) + (1 - a;^ cos $) {At^x^ cos $ - 2t^x'^) 
+ (cos $2 - x'^)t^ = 2t^x cos $ - t^x"^ + 2t^x. 

(27) 

It is easy to see from this equation that for the fully packed phase: tij oo, there 
are two possible choices of the background U(l) field's phase: cos$ = x"^ = ±1. In that 
limit eq. (P7D becomes 

(l_cos$) = i^. (28) 
In both cases we find that the spectrum 

cosh£(p) = 1 - 2(cos2p - 1)- + (cos4p + 4cos2p - 5) ^ (29) 

t t'^ 

when t — > oo becomes e{p) = for the all values of the momentum. 
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3 Representation of the partition function Z via two- 
and three-particle i?-matrices 

In this section we will demonstrate that the partition function (|^) can be constructed 
as a trace of the A^-th power (A^ is a size of the lattice in a time direction) of some transfer 
matrix T 

Z = trT^, (30) 

but unlike the case of ordinary integrable models, the transfer matrix T here is a product 
of four transfer matrices presented in the formulas (|^), (10) and ([TT|). However, if we 



look at the Fig. 3 under a 45° angle and use the corresponding time direction, we see 
only two different rows of products of constituent i?-matrices (see Fig. 4), corresponding 
to following Monodromy matrices in the so called braid formalism 



Uu) = rK-.i,3,(«)453,+1.3,+2(«), (31) 
j=0 
L 

Tl{u) = n^S!3,+l,3,+2(«)4?+2,3,+3(«)- (32) 
j=0 

Here the (m)'s are the ordinary two-particle scattering i?-matrices corresponding to 

the squares of the chain in the Fig. 4, and the -R-_i j j+i(u)'s correspond to the hexagons 
and represent three-particle scattering i?-matrices. Therefore the Monodromy matrix of 
the model is the product T{u) = Ti(m)To(m). 

The two-particle i?jj-matrix is an operator acting on the direct product of the two two 
dimensional spaces (spaces of the fermions with spin) defined on the sites i and j, which 



according to the technique developed in []T6| can be fermionised by considering the basis 
states as \k), k = 0,1, with |1) = c"'"|0): 

<^Ki)®|ji) = (^..)KK2)®|j2). (33) 



Graphically it is represented in Fig. 5b, in the same way as in [|1^] for the /^-matrices of 
the XXZ-model in the braid formalism. 

The three-particle i?jjfc-matrix is an endomorphism on the direct product of three 
two-dimensional spaces (see Fig. 5a) with a basis defined as before: 

® Ijl) ® l^l) = iR^Jk)'^^\^2) ® |J2) ® \k2). (34) 



Following to the [|T6| we can represent i^^^^-matrices (s = 2, 3) via fermionic creation- 



annihilation operators by considering the Hubbard operators 

X; = |j)(z|, I, J = 0,1, (35) 
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and taking into account that the fermionic Fock space is graded with the parities of the 
states defined as p{i) = i. Then, by definition 



and 

,(3) 



(36) 



R^L =R^Jk\^l)\Jl)\h){h\{J,\{^,\ 

= (i?.,.)SS?K2)|j2)|fc2)(fcl|0l|(.l| (37) 
— y ^) y^^jk)i-^j^ki^i2^j2^k2- 

The fermionic expression of for the XX Z model can be found in [|TB[(and refer- 
ences there), while the most general form of the three-particle scattering R^23 operator 
is 

^123 = Rooo^i^^ns + ^ool?^i^2^3 + ^ol[j^i^2% + i^lo^^inang + i^g^Jninang 

+R\'^\nin2n3 + R\l'^nin2n3 + R\l\nin2n3 + (^oo?4c3 + Roil4'^2)ni /oo\ 
+ {R\l'i4c3 + Rl%ctc2)n, + {Rh'olctcs + RTol4ci)n2 + (i?^l?c+C3+ ^''^ 
R'iw4ci)n2 + {R'oToctc2 + R%^4ci)n3 + (/?^?|c+C2 + R%\ctci)ns, 

Where = cfci and = 1 — rij. The two-particle i?|^'*-matrix of the XXZ model can 
be obtained from (^) by putting = 1 everywhere and taking {R'^'^^Yi^'jl = {R'^^^YiliyL 
and Rifjlo = otherwise. 

Both kinds of i?-matrices can be written as an exponent 

: expcfAijCj :, (39) 

(2) 

where the indices i,j run from 1 to 2 for R\2 operators and to 3 for the three-particle 
R^i23 operator. The notion : : means the normal ordering for the , Cj (i=l,3) and the 
hole ordering for the , C2 operators. 

The connection between the matrix elements Aij and the hopping parameters tj j in 
the action of the model can be found by the relation 

(^i|(V'2|(^3| : expc+AjjCj- : |V'3)|V'2)|V'i) = exp {-^piUj^pj). (40) 

and are 

4 _ J ^ 7^ i = 1,3 . _ j 1 + tjj i = 2 , s 

^^■"l 1 = 2 '^--1 -(l + t,,) z = l,3 • ^^'> 



In (40) the coherent states \ipi) are defined by the formulas (H). 
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The matrix elements (-R*-^-*)'^" in the expression ( PSD are connected with the Aij-s in 
f) by the following equations 



DlOl _ 1 DllO _ A plOl _ A 

-Kioi — -Kioi — ^23, -Kiio — ^32, 

/2[]l| = (l + An)(l-A22)+Ai2A21, i?^?l=v4i2, i??Jl=^21, 

R\l^, = (1 + A33)(l - A22) + A23A32, = Ai3(l - A22) + A12A23, 

= 1 - A22, R'lll = A3l(l - A22) + A21A32, 



(42) 



Koi = ^23(1 + ^11) - ^13^1, R'oih = ^32(1 + An) - A31A12, 

-^001 = ~^13 -RlOO = -^31, (43) 

= Ai2(l + A33) - A32A13, = A2i(l + A33) - A23A31, 

-^001 = 1 + ^11 ! 

= (1 - A22)(l + An + A33) - det A + A11A33 - A13A31 + A12A21 + A23A32, . . .^ 

DlOO — T _i_ A \^^) 
^100 — -L + ^33, 

RZ = {l + An){l + A,3)-Ar3A,,. 



(2) 

The two-particle i^^y -matrix can be obtained from these expressions by taking Ai^ = 
A^j = everywhere. 

(3) 

This is the general form of three-particle R-j^ matrix, but for the model under consid- 
eration d^) we should consider free fermionic limit and take A22 = 1 and ^13 = Ai-^ = 0. 

The condition of integrability of this model with two- and three-particle i?-matrices, 
namely the condition of commutativity of Transfer matrices ( pTD for different values of the 
spectral parameter, can be written as a modified Yang-Baxter equation in the following 
form 

Rl2{u,v)R234.iu)R4,5{u)Ri2{v)R234.{v) = -R234 (^^)-R45 (^^)-Rl2 (m)-R234 (m)R45 («, w) , (45) 

were Ilij{u,v) is the intertwiner operator. This equations definitely has a solution in free 
fermionic case corresponding to the model defined above, since we were able to diagonalise 
the Hamiltonian. It would be interesting to find a solution for a general case of R^^^ matrix 
(|M3) and i?(2)-matrix of the XX Z model. 
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Figure captions 

la. The BCC lattice 

lb. The Dual BCC lattice 

2. How the Manhattan Lattice is forming on surfaces on Dual BCC lattice 

3. Manhattan lattice appeared on 2D surfaces on Dual BCC lattice 

4. Two rows of Monodromy matrices 

5a. Graphical representation of - matrix 
5b. Graphical representation of matrix 
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Fig.l 




Fig.2 
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